Constraints on Cosmological Parameters from the Analysis of the Cosmic Lens All 
Sky Survey Radio-Selected Gravitational Lens Statistics 
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We derive constraints on cosmological parameters and the properties of the lensing galaxies from 
gravitational lens statistics based on the final Cosmic Lens All Sky Survey (CLASS) data. For a flat 
universe with a classical cosmological constant, we find that the present matter fraction of the critical 
density is f2 m = 0.311q'jI (68%) j^o xo (systematic). For a flat universe with a constant equation of 
state for dark energy w = p x (pressure) / p x (energy density), we find w < — 0.55^°'^ (68%). 

PACS numbers: 98.80.-k,98.80.Es,98.80.Cq,98.62.Py 



Gravitational lensing, in which a background source 
behind a foreground object is seen as distorted and 
(de)magnified image(s), is well-understood by the sim- 
ple physics of light deflections in a weak gravitational 
field. Hence, the analysis of lens statistics can provide 
a technique ]jj for constraining cosmological parameters 
that is both powerful and complementary to other meth- 
ods However, for lens statistics to be useful, it is vital 
to have an unbiased statistical sample of sources that is 
complete within well-defined observational selection cri- 
teria |j] and to properly take into account all factors in 
a statistical lensing model B . The statistical properties 
of gravitational lensing in a sample are the total lensing 
rate, the image separations, the lens redshifts, the source 
redshifts, and the image multiplicities. These properties 
depend not only on cosmological parameters but also on 
the properties of the lensing galaxies and the distribu- 
tions of the sources in redshift (z) and in luminosity in 
the observational selection waveband (e.g. radio here) j|. 

Recent advances in observations permit us to do a re- 
liable analysis of lens statistics to obtain limits on cos- 
mological parameters. First, the largest radio-selected 
galactic mass-scale gravitational lens search project to 
date, the Cosmic Lens All Sky Survey (CLASS), has now 
been completed, resulting in the largest sample suitable 
for statistical analyses Second, currently available 
independent data sets on the distributions of (flat spec- 
trum) radio sources in flux density and in redshift are 
concordant Q. Third, recent large-scale observations 
of galaxies, in particular the Two Degree Field Galaxy 
Redshift Survey (2dFGRS) and the Sloan Digital Sky 
Survey (SDSS), have produced converging results on the 



total galaxy luminosity function (LF; i.e., distribution of 
galaxy number density in luminosity), and there exist 
fairly reliable observational results that permit us to ex- 
tract from the total LF morphological type-specific LFs, 
i.e., an early-type (i.e. ellipticals and SO galaxies) LF and 
a late- type (i.e. spirals and others) LF j| In this Let- 
ter, we report the main results on cosmological param- 
eters from a likelihood analysis of lens statistics based 
on these crucial sets of data. The detailed procedure of 
the analysis will be described in a later publication Q 
in which extended analyses will be presented for broader 
range of parameters and with an emphasis on the global 
properties of galaxy populations. 

The CLASS statistical sample suitable for our analysis 
contains 8958 radio sources out of which 13 sources are 
multiply-imaged |3) . The lens search looked for multiple- 
lcnsing image components of the central compact radio 
core in the background sources. The CLASS statisti- 
cal sample is defined so as to meet the following ob- 
servational selection criteria j|: (1) the spectral index 
between 1.4 GHz and 5 GHz is flatter than —0.5, i.e. 
a > —0.5 with S v on v a where S u is flux density mea- 
sured in milli-jansky (mjy; 1 mjy = 10 -29 W m~ 2 Hz -1 ); 
(2) the total flux density of each source is > 30 mjy 
at 5 GHz; (3) the total flux density of each source is 
> 20 mjy at 8.4 GHz; (4) the image components in lens 
systems must have separations > 0.3 arcsec and the ratio 
of the flux densities of the fainter to the brighter com- 
ponent in double-image systems must be > 0.1. Number 
counts for sources at v = 5 GHz are consistent with a 
(broken) power-law differential number-flux density rela- 
tion \dN/dS\ (i.e. number of sources per flux density bin) 
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TABLE I: Gravitational lens systems in the well-defined 
CLASS statistical sample. z a and zi are source and lens red- 
shifts, respectively; for sources with unmeasured redshifts, we 
take z a = 2, which is the mean source redshift for the entire 
sample of CLASS lenses. The maximum image separation 
(A9) is given in arcsec. Image multiplicity is listed under 
Nj m . When the image splitting is due to multiple galaxies, 
the image separation and the image multiplicity are not used 
in our model fitting. The image separation of 2045+265 is 
not used because of uncertainties in the nature of the lensing 
scenario. Probable lens galaxy type identification (G-type) is 
coded as: e for an early-type and s for a spiral-type. 
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oc (S/S )- 7 ] with r] = 2.07 + 0.02 (1.97 + 0.14) for S > S° 
(S < S°) and S* = 30 mJy. These results are consis- 
tent with the prediction of the Dunlop & Peacock (1990) 
free-form model number 5 radio luminosity function that 
is consistent with redshift data at S > 1 mJy [|[ Q. 
Redshift measurements for a representative CLASS sub- 
sample and the prediction of the aforementioned Dunlop 
& Peacock (1990) model consistently suggest that the 
redshift distribution for the CLASS unlensed sources can 
be adequately described by a Gaussian model with mean 
redshift, (z) = 1.27 and dispersion 0.95 j| §. The prop- 
erties of the 13 multiply- imaged sources are summarized 
in Table 1. 

The luminosity function of galaxies is assumed to be 
well-described by a Schechter function H of the form 



dn 



exp(-L/L*), 



(1) 



where a is a faint-end slope and n» and L* are char- 
acteristic number density and characteristic luminosity, 
respectively. From Eq. (1) an integrated luminosity den- 
sity j is given by j = J °° dLL(dn/dL) — n a( L lr T{a + 2). 
The luminosity of a galaxy is correlated with its line- 
of-sight stellar velocity dispersion (cr) via the empirical 
relations, 



L 



^g0.4(M»-Af) 



(2) 



where M* and cr* are respectively characteristic ab- 
solute magnitude and characteristic velocity disper- 



sion, and 7 corresponds to the Faber- Jackson (Tully- 
Fisher) exponent for an early-type (late-type) popula- 
tion Q; throughout we take TFaber-jackson = 4.0 and 
TTuUy-Fisher = 2.9 §. The 2dFGRS team and the SDSS 
team independently determine the Schechter parameters 
for galaxies at z < 0.2 summed over all galactic types 
and their results are consistent with each other ||. The 
converging result can be expressed in the bj photometric 
system as M* — 51og 10 h — —19.69 ± 0.04 (hereafter h is 
the Hubble constant Ho in units of 100 km s _1 Mpc -1 ), 
a= -1.22+0.02, and n* = (1.71±0.07)x 10~ 2 h 3 Mpc" 3 , 
which gives an integrated luminosity density of j = 
(2.02 ± 0.15) x 10 8 h L Q Mpc' 3 . From the total LF 
we extract the type-specific LFs based on measured rel- 
ative type-specific LFs keeping the total luminosity den- 
sity fixed. The type-specific LFs are required because 
the early-type and the late-type populations are dynam- 
ically different and contribute to the multiple imaging 
in distinctively different ways, namely that the early- 
type population contributes more to the total lensing rate 
and on average generates larger image splittings Q . We 
use the relative LFs obtained from the Second South- 
ern Sky Redshift Survey (SSRS2), which is based on 
5404 local (z < 0.05) galaxies. We obtain the follow- 
ing early-type (e) and late-type (s) LFs expressed in the 
bj photometric system: for the early-type population, 
A/i e) - 5 log 10 h = -19.631^?, = -1.00 + 0.09, and 
ni e) = (0.64 + 0.19) x 10 -2 h 3 Mpc -3 ; for the late-type 
population, c*M = -1.22 + 0.02, Afi s) = -19.69 + 0.04, 



and 



,0) 



(1.20 ± 0.11) x 10" 2 h 3 Mpc- 3 . To break 



a degeneracy we have chosen = a and M* = Af* 
since the total LF is dominated by late-type galaxies. For 
alternative type-specific LFs obtained from the 2dFGRS, 
the reader is referred to j|, || . When the alternative type- 
specific LFs are used, our derived limits on cosmological 
parameters are essentially unchanged. 

The multiple imaging cross section of a galaxy can 
be related to and calculated from the line-of-sight ve- 
locity dispersion of the galaxy. However, the relation is 
not unique but depends on the mass profile and intrinsic 
shape of the galaxy. This dependence will be absorbed 
into a "dynamical normalization" factor A(/) below. It is 
known that projected surface densities of the inner cylin- 
drical regions of galaxies along the line-of-sight probed 
by gravitational lensing are well approximated by a sin- 
gular isothermal mass profile or a profile similar to it [Q . 
In this Letter, we adopt a singular isothermal ellipsoid 
(SIE) as a galaxy lens model, whose projected density 
can be expressed in units of lensing critical surface den- 
sity E cr = c 2 Dos/(47tGL'ol£ ) ls) as 



n{x,y) 



rcr VTKf) 



2 v^ + TV 



(3) 



where / is the apparent minor-to-major axis ratio, X(f) is 
the dynamical normalization factor for a given dynamical 
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model, and 



Finally, the integrated probability is given by 



47T 







DqlDls 
D os 



(4) 



is the critical radius for / = 1. £>ol, £>ls and Dos are 
angular diameter distances between observer (O), lens- 
ing galaxy (L), and source (S). Calculations of dynami- 
cal normalizations for axisymmetric cases (i.e. oblate and 
prolate) can be found in J|. For the case of equal num- 
bers of oblates and prolates, |A(/) — 1| < 0.06 for / > 0.4. 

Let s m (/) be the lensing cross section for image multi- 
plicity m (— 2, 4) due to the surface density given by 
equation (3) in units of r 2 r (Eq. 4). Then, we have 
s m(f) = [A(/)] 2 s m (/) where s m (/) can be evaluated nu- 
merically Q. The differential probability for a source 
with redshift z s and flux density 5* to be multiply-imaged 
with image multiplicity m and image separation A9 to 
A9 + d(A9), due to a distribution of intervening galaxies 
of type g (g = e, s) at z to z + dz modeled by SIEs and 
described by a type-specific Schechter LF, is given by 
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V c) AOA Ad, 



DqlDls 
Dos 

(q 7 +7+2)/2 



x exp 
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7/2' 



B m {z„S), (5) 



where B m (z s ,S) is a magnification bias factor (i.e., an 
over-representation of multiply-imaged sources because 
of flux amplifications) which depends on a magnification 
probability distribution and \dN/dS\ and can be calcu- 
lated numerically [Q. In Eq. (5), AO, is a characteristic 
image (angular) separation given by 
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w&(t)' 



Dos 
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The differential probability with A9 to A9 + d(A9) can 
be obtained from Eq. (5): 



dp 



(s) 
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dp 



d{A9) 7 dzd{A9) 



dz. 



The differential probability with any A9 > A9 
galaxies at z to z + dz is given by 



(7) 

lilin due to 



dp 
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dz 



s m (f)Ww'T l a + l + - 
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dt 



xn. 



DolD ls 



cJ \ Dos 

A9 min d 2 p (9) 
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For the sample of CLASS sources containing Al 
multiply-imaged sources and N\j unlensed sources, the 
likelihood £ of observing the number and the properties 
of the multiply- imaged sources (Table 1) and the num- 
ber of the unlensed sources given the statistical lensing 
model, is defined by 



N v 

ln£ = ^ln 



1 - 



m=2A 



k=l 



(10) 

where the integrated probability (Eq. 9) Pm l \i) = 
Pm\i) + Pm\i) and the differential probability (Eqs. 5, 
7, or 8) Spm ne \k) is given by 

8pt nc) (AO = w (e) (k)5p$ (k) + (k)6p£ (fc) (11) 

with w^(k) + w^(k) = 1. If the A:-th lensing galaxy 
type is known to be an early- type (late- type), (k) = 1 
(w^ (k) = 1). If the fc-th lensing galaxy type is unknown, 
we use w^)(k) = 5p${k)/[8p { $(k) + 6p$(k)] if the lens 
redshift is known, and we use w^ e ' (k) = 0.8 and w^ s ' (k) — 
0.2 otherwise. From Eq. (10), a "x 2 " is defined as 



X 



-21n£. 



(12) 



We calculate the x 2 (Eq. 12) over grids of the model 
free parameters (see below; in each grid point the \ 2 is 
minimized over the nuisance parameters) and determine 
the best-fit values and errors of the free parameters. We 
allow most of the critical parameters to be free while 
we fix the parameters that are well-determined by ob- 
servations or to which the x 2 is insensitive. The free 
parameters to be determined by minimizing the x are 
as follows: present matter density Q m , present vacuum 
density Oa or dark energy density Q x (all expressed as 
fractions of the present critical density) with its constant 
equation of state w = p x j p x (where p x and p x are its 
isotropic pressure and uniform energy density), early- 
type and late- type characteristic velocity dispersions ai^ 

(s) 

and cr* , and the mean apparent axial ratio of galaxies 
/. A critical but fixed parameter is the early-type char- 
acteristic number density n^ . However, whichever of 
the presently available choices of type-specific LFs we 
use, the derived constraints on cosmological parameters 
turn out to be virtually the same. One key assumption 
we make is that the (relatively) locally determined early- 
type number density nt^ along with the faint-end slope 
have not evolved since z ~ 1. This appears to be a 
valid assumption based on current evidence and under- 
standing j7j. However, if there is any significant reduction 
in the early-type number density at z ~ 0.6 compared 
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with the present epoch, our presently determined mat- 
ter density (dark energy density) will be overestimated 
(underestimated) . 

We present our results based on the type-specific LFs 
derived using the SSRS2 LFs and a dynamical normaliza- 
tion assuming equal frequencies of oblates and prolates. 
Figure 1 shows likelihood regions in a matter plus vac- 
uum cosmology: we find the combination Q\ — 1.2f2 m = 
0.40±g;g| (68%) (±5:lf at 95%). Figure 2 shows the likeli- 
hood in the w-fl m plane with the prior that the universe 
is flat; we find w < -0.55±^ (68%). This limit is in 
agreement with other recent results ||. For a flat uni- 
verse with a cosmological constant (i.e. for w = —1), we 



find fL 



0.31 



+0.27 
-0.14 



(68%) (+0.2I at 95%) 



+0.12 
-0.10 



(sys- 



tematic). The identified systematics are due to the un- 
certainties in the redshift distribution and the differen- 
tial number-flux density relation of the sources. For the 
flat w = — 1 case, other fitted values of the parameters 
are: a[ e) = 198±p km s"\ a[ s) = 117±|f km s" 1 , and 
/ < 0.83, all at 95%. 

In conclusion, the likelihood function (or, x 2 ; Eq. 12) 
of the CLASS data has well-defined confidence regions 
in the parameter space defined by our cosmological and 
galactic parameters. Our constraint in the f2 m -fiA plane 
based solely on lens statistics agrees with that from 
Type la supernovae observations Our determined 
value of f2 m in flat cosmologies is also in good agree- 
ment with several recent independent determinations H . 
Compared with previous works in lens statistics our 
results imply not only positive but higher dark energy 
density. Therefore, our results, which are based on dif- 
ferent physics, assumptions and astronomical relations, 
and independent sets of data, add to the evidence for a 
universe with low matter density and high dark-energy 
density. 
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